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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 
RESEARCH  MEMORANDUM 

MEASUREMENTS  OF  NORMAL-FORCE-COEFFICIENT  FLUCTUATION  ON 
FOUR  9-PERCENT -THICK  AIRFOILS  HAVING  DIFFERENT 
LOCATIONS  OF  MAXIMUM  THICKNESS 
By  Milton  D.  Humphreys 

SUMMARY 


A two-dimensional  wind-tunnel  investigation  of  the  effects  of 
maximum-thickness  location  on  fluctuating  pressures  and  normal-force 
coefficients  on  9-percent -thick  airfoils  indicated  that  for  normal-force 
coefficients  up  to  0.6  the  section  variables  had  little  effect  on  the 
pressure  pulsations  and  root -mean- square  normal-force-coefficient  fluc- 
tuations. At  normal-force  coefficients  above  0.6,  up  to  approximately 
the  maximum  normal-force  coefficient,  considerable  reductions  in  root- 
mean-square  normal-force  fluctuations  were  obtained  by  using  an  airfoil 
having  the  maximum  thickness  located  at  or  behind  the  42. 5-percent -chord 
station.  From  the  standpoint  of  local  panel  loads,  none  of  the  airfoils 
had  consistently  lower  fluctuating  loads  throughout  the  speed  and  normal- 
force-coefficient  ranges.  It  appeared,  however,  that  the  NACA  65AOO9  air- 
foil would  probably  be  the  best  choice  of  the  sections  tested  for  oper- 
ation throughout  the  high  subsonic  Mach  number  range. 


INTRODUCTION 


Former  investigations  (for  example,  refs.  1 and  2)  have  provided 
useful  information  on  the  character  and  magnitude  of  the  local  pressure 
pulsations  on  airfoils  at  high  subsonic  speeds.  These  results  can  be 
used  for  the  prediction  of  local  panel  loads,  but  the  fluctuating  normal- 
force  coefficient  of  the  section,  which  is  of  more  significance  in  buf- 
feting studies,  cannot  be  obtained  frcm  the  pressure-pulsation  records 
unless  detailed  consideration  is  given  to  the  phase  relationships  of  the 
pressures  at  different  chordwise  stations.  In  the  present  study  use  is 
made  of  an  electrical  integrating  device  to  obtain,  from  the  pulsating- 
pressure  measurements,  the  desired  root -mean- square  values  of  the  fluc- 
tuating normal-force  coefficient  for  several  9-percent -thick  airfoils. 


CONFIDENTIAL 


2 


CONFIDENTIAL 


NACA  PM  L54B22 


Previous  pres sure -pulsation  studies  for  6-percent -thick  sections 
(ref.  1)  indicated  considerably  lower  pulsations  for  the  NACA  65A006  sec- 
tion than  for  the  NACA  64A006  section.  In  the  present  invest igation, 
therefore,  a range  of  maximum-thickness  positions  was  covered  by  investi- 
gation of  the  NACA  63AOO9,  64A009,  65AOO9,  and  16-OO9  sections.  The  tests 
covered  a Mach  number  range  of  0.5  to  1.0,  angles  of  attack  of  0°  to  12°, 
and  a Reynolds  number  range  of  1.5  x 106  to  2.1  x 106.  The  results  of  a 
somewhat  similar  investigation  of  different  profiles  are  reported  in 
reference  3. 


SYMBOL 


M free-stream  Mach  number 

q free-stream  dynamic  pressure,  lb/sq  ft 

Ap  double-amplitude  pressure  pulsation,  lb/sq  ft 

cn  section  normal-force  coefficient 

ca  section  drag  coefficient 

375  chordwise  average  of  Ap/q 

Acn  square  root  of  mean  square  normal-force-coeff icient 

fluctuations  (unsteady  component  of  the  normal-force 
coefficient) 

a angle  of  attack,  deg 


APPARATUS,  MODELS,  AND  TESTS 


The  investigation  was  conducted  in  the  Langley  4-  by  19-inch  semiopen 
tunnel  (fig.  1).  Dried  air  at  20  pounds  per  square  inch  absolute  flowed 
from  the  settling  chamber  through  the  test  chamber  and  was  discharged  to 
the  atmosphere  through  the  diffuser.  The  test-section  Mach  number  was 
regulated  by  a variable-area  throat  located  in  the  diffuser  downstream 
from  the  test  section.  The  tunnel  calibration  and  operation  were  similar 
to  that  described  in  reference  4. 

Jet-boundary  corrections  for  this  type  of  test  have  not  yet  been 
determined  at  high  subsonic  Mach  numbers;  therefore  no  correction  has 
been  applied  to  any  of  these  data.  For  incompressible  potential  flow 
the  correction  to  the  angle  of  attack  is  the  major  correction  and  is 
given  (for  the  tunnel  configuration  used)  by  ■ a^est  - 1.85cn. 
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Jet-boundary  effects  are  discussed  in  more  detail  in  references  4 and  5, 
where  it  is  indicated  that,  except  for  the  angle-of -attack  correction, 
the  jet-boundary  effects  are  probably  not  large. 

The  4-inch-chord  models  investigated  and  their  position  of  maximum 
thickness  are  as  follows: 


Airfoil 

Position  of  maximum  thickness. 

percent  chord 

NACA  65A009 

36.25 

NACA  64A009 

59-0 

NACA  65A009 

42.50 

NACA  16-009 

50.0 

The  model  ordinates  are  given  in  reference  4.  Pressure  orifices  were 
located  on  the  upper  surface  of  the  models  at  the  5.1-,  14-,  25-,  37*5-> 
50- , 62.5-,  75-,  and  87 .5-percent -chord  stations.  The  lower-surface  flow 
conditions  were  not  investigated,  since  an  earlier  investigation  (ref.  2) 
showed  that  under  lifting  conditions  the  major  flow  unsteadiness  was  on 
the  upper  surface  of  the  airfoil.  This  may  lead  to  inaccuracies  in  the 
force  fluctuations  at  angles  of  attack  near  zero  liftj  however,  in  this 
region  the  fluctuations  are  small  and  inaccuracies  of  little  significance. 

Miniature  electrical  pressure  gages  and  a recording  system  having  a 
flat  frequency  response  from  40  to  500  cycles  per  second  were  used  to 
measure  the  fluctuating  pressures  at  each  of  the  eight  pressure  orifices 
along  the  chord  of  the  airfoils.  In  the  first  phase  of  the  investigation 
the  electrical  output  of  these  gages  was  recorded  to  give  the  time  history 
of  the  fluctuating  pressures  (as  in  fig.  2).  From  these  records  a typical 
or  predominant  double-amplitude  value  of  the  fluctuating  pressure  Ap  was 
chosen  by  inspection  (as  in  ref.  1)  and  plotted  to  give  the  distribution 
of  Ap/q  along  the  airfoil  chord  (fig.  5).  The  data  thus  obtained  are 
useful  in  the  determination  of  predominant  fluctuations  of  local  panel 
loads . 

In  the  second  phase  of  the  investigation  the  electrical  outputs  of 
the  same  miniature  electrical  pressure  gages  were  added  simultaneously  by 
means  of  an  electrical  integrator,  a dynamic  version  of  the  static  inte- 
grator described  in  reference  6,  to  produce  an  electrical  signal  corre- 
sponding to  the  instantaneous  fluctuating  component  of  the  normal  force. 
These  signals,  over  a period  of  time,  were  converted  by  a vacuum  thermo- 
couple and  millianmeter  into  a meter  reading  corresponding  to  the  square 
of  the  mean  fluctuating  normal  force.  This  information  was  easily  trans- 
posed into  the  root  mean  square  of  the  fluctuating  normal  force  (fig.  4). 

A laboratory  calibration  of  the  system  indicated  that  the  response  was 
approximately  1 to  1 up  to  500  cycles  per  second  and  was  amplified  con- 
siderably (up  to  1.8  times  at  atmospheric  pressure)  between  500 
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and  1,000  cycles  per  second.  The  amplification  was  reduced  considerably 
as  the  applied  static  pressure  was  reduced  with  increasing  Mach  number; 
hence,  the  component  of  the  data  contributed  by  the  higher  frequencies 
became  more  accurate  at  the  higher  Mach  numbers.  No  correction  for  fre- 
quency response  has  been  made  to  these  data,  since  laboratory  investiga- 
tions have  shown  that  only  negligible  phase  lag  existed  in  the  measuring 
apparatus.  Further,  since  all  the  eight  gages,  all  orifices,  and  all 
tubing  were  practically  identical,  their  phase  characteristics  would  be 
identical  and  would  have  no  effect  on  the  integrated  results. 

To  check  the  reproducibility  of  the  root -mean-square  normal-force 
fluctuation  for  a given  test  condition,  three  separate  runs  were  made  — 
two  using  the  electrical  integrator,  and  a third  run  from  which  the  pres- 
sure record  of  figure  2 was  mechanically  integrated.  The  results  of  the 
mechanical  integration  yielded  a root-mean-square  normal-force-coefficient 
fluctuation  of  0.055,  which  compares  favorably  with  electrically  inte- 
grated values  from  two  other  tests  of  0.043  and  0.060  for  this  test  point 
(fig.  5) • The  repeatability  for  three  separate  tests  obtained  at  a 

Reynolds  number  of  1.9  x 10^  and  at  maximum  normal-force  coefficient  is 
excellent,  in  view  of  the  very  unsteady  nature  of  the  flow  phenomenon 
measured  at  this  test  condition. 

Pressure-pulsation  data,  root -mean-square  force -fluctuation  data, 
and  schlieren  motion  pictures  were  obtained  at  angles  of  attack  from  0° 
to  12°.  The  Mach  number  ranged  from  0.5  to  1.0  and  the  corresponding 
Reynolds  numbers  were  approximately  1.5  x 106  to  2.1  x 106. 


RESULTS  AND  DISCUSSION 
Pressure  Fluctuations 


Pressure-pulsation  records,  as  shown  in  figure  2,  were  used  to  obtain 
typical  or  predominant  values  of  the  double -amplitude  pulsating  pressures 
expressed  as  Ap/q.  The  chordwise  variation  of  pressure  pulsations  on 
the  9-percent -thick  airfoils  at  several  constant  normal-force  coefficients 
is  presented  in  figure  3 to  show  the  effect  of  maximum-thickness  location 
and  Mach  number  on  the  pressure  pulsations.  At  low  normal-force  coeffi- 
cients of  0.2  and  0.4,  the  maximum  pressure  pulsations  are  less  than 
10  percent  of  the  dynamic  pressure  and  the  effect  of  the  location  of  maxi- 
mum thickness  on  pressure  pulsations  is  negligible.  At  normal-force  coef- 
ficients of  0.6  and  0.7  the  maximum  pressure  pulsations  become  significant, 
reaching  values  of  approximately  30  percent  of  the  stream  dynamic  pressure. 
While  the  NACA  63A009  airfoil  has  the  lowest  values  of  Ap/q  at  a Mach 
number  of  0.5,  it  shows  progressively  higher  values  as  the  Mach  number  is 
increased  and  has  the  highest  level  of  pressure  pulsations  at  Mach  numbers 
of  0.9  and  1.0.  The  NACA  64A009  airfoil  has  moderately  high  levels  of 
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Ap/q  at  all  Mach  numbers , without  ever  having  the  lowest  values.  The 
NACA  65AOO9  airfoil,  while  having  high  pulsation  levels  at  a Mach  number 
of  0.5>  shows  the  lowest  level  at  Mach  numbers  from  0.6  to  0.8  and  has 
next  to  the  lowest  levels  above  a Mach  number  of  0.8.  The  NACA  16-009  air- 
foil exhibits  high-level  pulsations  at  Mach  numbers  from  0.5  to  0.8,  but 
has  the  lowest  level  at  0.9  and  1.0.  Although  none  of  the  airfoils 
investigated  had  consistently  lower  pressure-pulsation  levels  at  all  Mach 
numbers,  it  appears  that  the  NACA  65A009  airfoil  would  probably  offer  the 
best  compromise. 


Normal-Force-Coefficient  Fluctuations 

Figure  2 shows  a sample  record  of  the  pressure  pulsations  on  an 
NACA  64A009  airfoil  at  10°  angle  of  attack  and  a Mach  number  of  0.77  and 
illustrates  the  phase  differences  that  can  occur.  Phase  differences  of 
as  much  as  180°  between  the  forward  gage  and  the  rear  gage  are  apparent 
in  figure  2,  and  these  large  phase  differences  are  responsible  for  the 
fact  that  the  root -mean-square  fluctuating  normal-force  coefficient  Acn 

is  of  smaller  magnitude  than  estimates  of  this  coefficient  made  from  the 
pressure  data  without  regard  for  phase  relations.  Further,  these  phase 
differences  may  produce  fluctuations  in  pitching-moment  coefficient. 

Before  examining  the  data  for  normal -force  fluctuations  in  more 
detail,  it  is  of  interest  to  compare  the  general  level  of  the  force  fluc- 
tuations with  the  level  of  the  predominant  pressure  fluctuations.  In 
order  to  make  such  a comparison  the  Ap/q  diagrams  of  figure  5 were 

integrated  to  obtain  an  average  value  Ap/q  along  the  chord.  This 
chordwise  average  value,  of  course,  is  arrived  at  without  consideration 
of  phase  relations  or  wave  form.  Furthermore  it  is  a double-amplitude 
value  by  definition.  These  factors  all  tend  to  make  Ap/q  several  times 
larger  than  the  measured  root -mean- square  normal-force -coefficient  fluc- 
tuation Acn. 

Figure  4 shows  a comparison  of  Ap/q  with  the  root -mean-square 
fluctuating  normal-force  coefficient  Acn  for  the  NACA  64A009  airfoil 
at  several  angles  of  attack  over  a range  of  Mach  numbers.  The  curves 
for  ^p/ q have  values  approximately  5 times  those  of  the  curves  for  Acn 
for  angles  of  attack  up  to  6°  over  most  of  the  Mach  number  range.  At  8° 
and  10°,  the  values  of  the  Ap/q  curves  increase  to  7 times  those  for 
the  Acn  curves.  The  increase  in  this  ratio  from  5 to  7 when  going  from 

low  to  high  angle -of -attack  ranges  is  primarily  attributed  to  an  increase 
in  the  phase  difference  of  the  pressures  along  the  chord,  since,  at  the 
low  angles  of  attack  (up  to  6°),  the  pressures  have  been  observed  to  be 
more  closely  in  phase  (ref.  2)  than  at  angles  above  6°,  where  phase 
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differences  of  as  much  as  l80°  have  frequently  been  observed  (for  example, 
fig.  2). 

As  previously  explained,  the  Ap/q  curves  in  figure  4 show  double- 
amplitude values  and  thus  these  data  should  be  reduced  by  a factor  of  2 
for  direct  comparison  with  the  curves  for  the  root -mean-square  fluctu- 
ating normal-force  coefficient.  When  this  is  done,  the  Ap/q  data  are 
still  approximately  3 times  the  magnitude  of  the  Acn  data.  This  resid- 
ual factor  includes  effects  of  phase  differences,  wave  form,  and  also  an 
unknown  factor  dependent  upon  the  arbitrary  method  used  to  select  the 
predominant  Ap/q  values.  It  is  notable  that  the  general  level  of  the 
normal-force-coefficient  fluctuations  is  low  — even  at  high  angles  of 
attack  and  at  Mach  numbers  between  0.7  and  0.9  the  normal -force  fluctu- 
ations seldom  exceed  5 percent  of  the  mean  normal  force. 

The  variation  of  the  root -mean-square  fluctuating  normal-force  coef- 
ficient Acn  with  Mach  number  is  presented  in  figure  5*  The  general 
trends  of  Acn  for  the  airfoils  are  in  agreement  with  the  Ap/q  data 
shown  in  figure  3.  For  angles  of  attack  up  to  6°  and  for  all  Mach  num- 
bers, the  variations  of  Acn  with  Mach  number  are  approximately  equal 
for  all  the  airfoils,  with  the  NACA  63A009  and  65A009  airfoils  generally 
having  slightly  lower  values  than  the  other  two  airfoils.  Increasing 

the  angle  of  attack  to  8°  produces  only  a negligible  change  in  the 
Acn  levels  for  Mach  numbers  up  to  0.725  and  this  is  attributed  to  the 
steady  nature  of  the  separated  flow  observed  on  all  the  airfoils 
(fig.  6(a)).  Above  a Mach  number  of  0.725,  with  the  flow  attempting  to 
attach  at  the  nose,  the  6-series  airfoils  show  rapid  rises  in  the  normal- 
force  -coefficient  fluctuations.  The  flow  attachment  occurs  on  the  NACA 
l6-series  airfoil  without  much  unsteadiness  and  a relatively  steady 
normal-shock  formation  is  obtained,  with  a steady  separated  flow  aft  of 
the  shock,  resulting  in  the  generally  low  root -mean- square  normal-force- 
coefficient  fluctuations  shown  for  the  NACA  l6-series  airfoil.  At  angles 
of  attack  above  8°,  the  behavior  of  the  curves  for  all  the  airfoils  is 

similar  to  that  at  6°,  except  that  the  levels  of  Acn  are  much  higher 
than  at  lower  angle s of  attack  for  all  Mach  numbers,  and  the  Acn  rise 
is  more  abrupt.  The  Acn  levels  do  not  consistently  rise  or  fall  with 
a systematic  variation  of  the  position  of  maximum  thickness,  the  NACA 
64A009  airfoil  generally  having  the  highest  value  of  Acn  over  the  Mach 
number  range.  A large  rise  in  Acn  occurs  at  angles  of  attack  of  10° 
and  12°  at  Mach  numbers  around  0.7  and  is  produced  by  the  oscillations 
of  shock  and  separation  point  on  the  airfoils  accompanying  the  unsteady- 
flow  attachment  (fig.  6(b)).  Lower  root -mean- square  normal-force- 
coefficient  fluctuations  are  shown  for  the  NACA  65A009  and  16-009  air- 
foils because  of  milder  shock  and  separation-point  oscillations  and  less 
violent  flow  movement  than  for  the  other  airfoils  at  10°  angle  of  attack. 
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Figure  7 shows  the  variation  of  the  Mach  number  for  the  Acn  rise 

with  angle  of  attack  and  the  Mach  number  for  the  normal-force  break  and 
drag  rise.  There  is  no  absolute  correlation  between  these  three  parame- 
ters because  of  the  erratic  variation  of  the  Mach  number  for  the 
Acn  rise  with  regard  to  the  other  parameters.  The  Acn  rise,  however, 
generally  occurs  at  higher  Mach  numbers  than  those  for  the  drag  rise  but 
at  a Mach  number  close  to  that  for  the  lift  break. 

The  variations  in  Acn  with  cn  for  the  four  airfoils  at  various 
Mach  numbers  are  shown  in  figure  8.  At  Mach  numbers  of  0.5  to  0.85,  Acn 
is  not  greatly  affected  by  changes  in  normal-force  coefficient  up  to 
cn  = 0.6,  but  above  cn  = 0.6,  at  values  approaching  the  maximum  lift 
coefficient  for  the  sections,  there  is  a rapid  rise  in  Acn.  Similar 

effects  were  noted  in  the  variation  of  Ap/q  with  Mach  number  and  normal- 
force  coefficient  (fig.  3)*  When  the  Mach  number  is  increased  above  0.85, 
the  rapid  rise  in  Acn  at  the  higher  lift  coefficients  gradually  dimin- 
ishes and  at  Mach  numbers  of  the  order  of  0.98  there  is  little  change  in 
Acn  with  increasing  normal-force  coefficient  throughout  the  range.  These 
results  indicate  that  in  the  region  of  small  changes  in  Acn,  which  corre- 
sponds generally  to  normal-force  coefficients  below  0.6,  there  is  little 
effect  of  section  on  the  level  of  force  pulsation  except  at  a Mach  num- 
ber 0.5,  where  the  NACA  63A009  and  65A009  airfoils  have  lower  levels.  In 
this  region  of  normal-force  coefficients  approaching  maximum  lift,  where 
the  fluctuating  load  is  increasing  rapidly,  the  values  of  Acn  measured 
on  the  NACA  65AOO9  airfoil  were  smaller  than  those  measured  on  the  other 
airfoils.  Similar  effects  are  noted  in  figure  9,  where  the  normal-force 
coefficient  and  Mach  number  at  which  the  pulsating  normal  force  exceeds 
the  specific  values  of  0.012,  0.014,  and  0.016  are  shown  for  the  various 
sections.  These  results  also  show  no  consistent  effect  of  changes  in 
maximum-thickness  location;  however,  at  high  Mach  numbers  the  NACA 
65AOO9  airfoil,  in  general,  showed  high  cn  levels  for  specific  values 
of  Acn  throughout  the  Mach  number  range. 

Figure  10  is  a summary  plot  presenting  the  effects  of  maximum- 
thickness location  for  three  6-series  airfoils  and  one  unrelated  l6-series 

airfoil  on  Acn  for  several  constant  normal -force  coefficients.  The  var- 
iation of  Acn  with  maximum-thickness  location  is  characterized  by  peak 

values  of  Acn  occurring  around  the  39-percent-chord  maximum-thickness 

location  (NACA  64A009  airfoil),  while  the  airfoils  using  a forward  or  a 
more  rearward  maximum-thickness  location  generally  show  lower  values  of 

Acn.  The  airfoils  having  a maximum-thickness  location  at  or  behind  the 
k-2. 5-percent -chord  station  (NACA  65A009  and  16-009)  generally  had  the 

lowest  levels  of  Acn,  except  at  cn  = 0.2,  where  the  levels  of  Acn  are 
of  limited  significance. 
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CONCLUDING  REMARKS 


A two-dimensional  wind-tunnel  investigation  of  the  effects  of  the 
location  of  maximum  thickness  on  pressure  pulsations  and  root -mean- square 
force  fluctuations  on  three  6-series  airfoils  and  one  l6-series  airfoil 
of  9-percent  thickness  indicated  that: 

1.  Marked  variations  were  found  in  the  phase  relationships  of  the 
pressures  at  the  various  chordwise  stations.  Thus,  the  average  pulsating 
pressure  along  the  chord  cannot  be  used  as  an  average  value  of  the  fluc- 
tuating normal  force.  Frcm  the  present  results  it  appears  that  the  cor- 
rect root -mean- square  value  of  the  normal-force-coefficient  fluctuations 
is  of  the  order  of  one-third  of  the  value  estimated  from  a chordwise 
average  of  the  predominant  pressure  fluctuations  by  assuming  all  of  the 
pressure  fluctuations  to  be  in  phase.  The  results  also  indicated  that 
phase  differences  tend  to  become  more  pronounced  at  the  higher  angles  of 
attack. 

2.  The  largest  normal-force  fluctuations,  which  occurred  in  the  Mach 
number  range  from  0.7  to  0.9  at  the  highest  test  angles  of  attack,  seldom 
exceeded  about  5 percent  of  the  mean  normal  force.  A comparison  of  the 
root -mean- square  force -fluctuation  levels  for  four  9-percent-thick  air- 
foils (NACA  63AOO9,  64A009,  65AOO9,  and  16-009)  revealed  no  marked  dif- 
ferences at  normal-force  coefficients  up  to  about  0.6.  At  higher  values 
of  normal-force  coefficient,  however,  the  NACA  65AOO9  and  NACA  16-009  air- 
foils, which  had  more  rearward  maximum-thickness  positions  than  the  other 
two  sections,  had  lower  levels  of  fluctuating  normal  force. 

3.  From  the  standpoint  of  local  peine  1 loads,  none  of  the  airfoils 
had  consistently  lower  loads  throughout  the  speed  and  normal -force - 
coefficient  ranges.  It  appeared,  however,  that  the  NACA  65AOO9  airfoil 
would  probably  be  the  best  choice  of  the  sections  tested  for  operation 
throughout  the  high  subsonic  Mach  number  range. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va .,  February  10,  195^ • 
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Expanded  section  of  record 


Figure  2.-  Topical  pressure  pulsations  on  an  NACA  64A009  airfoil,  illus- 
trating phase  differences  that  occur,  a = 10°;  M = O.77. 
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(a)  cn  = 0.2.  (b)  cn  = 0.4. 

Figure  5 •-  The  effect  of  maximum- thickness  location  and  Mach  number  on 
the  chordwise  pressure  pulsations  on  9-percent-thick  NACA  airfoils 
at  several  constant  normal-force  coefficients. 
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(c)  cn  = 0.6.  (d.)  cn  = 0.7. 


Figure  3*-  Concluded. 
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(a)  a = 10°;  M = O.65. 

Figure  6.-  Schlieren  motion-picture  photographs;  250  frames  per  second.. 
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L- 83292.1 

NACA  67A009  NACA  65AOO9 

(b)  a = 10°;  M = 0.77'. 

Figure  6.-  Concluded. 
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